Dissipation-induced coherent structures in Bose-Einstein condensates 
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Introduction. - Dissipation is one of the main forces act- 
ing against the formation of nonlinear coherent structures 
in extended systems. When dissipation is present in sys- 
tems without additional gain mechanisms, typically all 
excitations decay into the regime of linear waves. In 
this Letter we discuss how, contradicting this general 
principle, a localized dissipation can be used to engi- 
neer the phase of systems governed by complex order pa- 
rameters through the generation of currents that imprint 
the required phases in the system. Although our results 
have general implications we discuss examples of systems 
ruled by an universal model of mathematical physics: the 
nonlinear Schrodinger equation (NLS). As to application 
fields we will focus on studying the possibility of control- 
ling the phase of Bose-Einstein condensates (BECs) to 
generate different types of coherent structures. 

It was soon after the realization of BECs that their 
potential to support nonlinear coherent structures was 
recognized. The list of experimentally found excitations 
includes: dark P, 0] and bright [3, [3 solitons, vortices, 
vortex lattices and related structures[l| , vortex rings 0] , 
gap solitons 0, shock waves different types of vec- 
tor solitons and Faraday waves [13] . Many techniques 
have been discussed to generate these structures but es- 
sentially all of them are "conservative" in nature based 
either on time- varying potentials, spatially selective op- 
tical transitions or on tailored interatomic interactions. 

In ultracold quantum gases dissipative mechanisms are 
related to inelastic collisions 
thermal component (TJ" 
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interaction with the 
or collapse dynamics Q, \\^ . 
While the latter can result in surviving matter-waves 
leading to a "non-destructive" effect of dissipation 
on certain nonlinear structures, dissipation is generally 
found to damp the excitations and act against the gen- 
eration or survival of coherent structures. In this Letter, 
however, we show that under certain conditions a prop- 
erly localized dissipation, e.g. the one recently demon- 
strated with the help of a focused electron beam 
can he used for generation and control of matter waves. 



Physical system and model equations.- We consider a 
quasi-one dimensional (ID) BEC in the mean field limit 
subject to a localized dissipative defect JJ| described 
phenomenologically by the distribution ^{x): 
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and 
are the 



Here time and space are measured in units of 2 /to 
a± = y''h/muj± respectively, where ll>j_ and a_ 
frequency and linear oscillator length of the transverse 
harmonic trap, and a =sign(as). The macroscopic wave 
function is normalized as J \ip\'^dx = 4J^o\as\/a±, A/q 
and fls being the initial number of atoms and the scat- 
tering length, respectively (see e.g. [lB|)- The applica- 
bility of the ID limit is determined by the condition 
a± <^ — [Snnolas])^^^^ ^ a||, where a|| is the longitu- 
dinal dimension of the condensate, ^ is either the healing 
length, when excitations against a constant background 
of density ng are under consideration, or a soliton width 
when one deals with a spatially localized solution (then 
nn is a maximal density of atoms). The local dissipation 
[iSj l is characterized by a temporal scale St which is the 
time interval between subsequent events of eliminating 
individual atoms from the atomic ensemble. The mean- 
field approximation is applicable if St is negligible, i.e. 
if uj±St <^ 1, what is verified for typical configurations 
where u!± ^ 2tt x lOOHz and St ~ 100 /iS. 

Effect of the localized dissipative defect. - To under- 
stand the physics of the phenomena discussed below we 
notice that the rate of particle decay in ([T]) is given by 
j{x)n{x) where n is the local particle density. Thus, the 
effect dissipation does not act close to the zeros of n{x) 
and eventually the loss can become negligible on the scale 
of the total population of the ground state provided the 
location of the defect matches the zeroes of the density. 

Let us first concentrate on repulsive interatomic inter- 
actions (tJ = 1). If we define C to be the defect size, i.e. 
the domain where j{x) significantly deviates from zero 
it is obvious that when C ^ then j{x) — joS{x), S{x) 
being the Dirac delta, and the real constant 70 character- 
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izing the strength of the defect (a dissipative (5- impurity). 
In that hmit a stationary dark sohton solution of Eq. ([T]) 
has the form ip^i = pexp (— zp^t) tanh (yoa;/-\/2). Obvi- 
ously, in that limit the dissipative losses, being localized 
exactly on the zeros of the density, are negligible. 

In the opposite limit where 7(0;) = 70 =const a simple 
solution of Eq. ^ reads 

^{x,t) = poe'''"'-^'^'\ ^ |L (1 „ e-^-'o*) (2) 

where po is the initial amplitude. 

Let us now consider the more interesting and physically 
relevant case of a finite size dissipative "defect" 

7(a:) =7oexp(-a:V2C'). (3) 

The decay of the atomic density in the center of the de- 
fect is characterized by the time = I/70. The sec- 
ond temporal scale depends on the size of the defect 
and is given by the time tg required for sound waves 
to leave the defect. Since the speed of sound near the 
defect core is c{t) = V^p{t) « y/2poe^'^°* we estimate 
ts = - log [1 - 7oC/(2\/2po)] ■ Thus, if U < t^, sound 
waves can emerge from the defect [Fig. (Tfc)]. However if 
tg 3> t-y the sound waves are strongly attenuated and the 
density decays smoothly in the vicinity of the defect. As 
it is clear this last scenario is always realized for C ^ Cc 
where Cc = 2V2po/7o [Fig. [Hd)] . 

To confirm these ideas, we have simulated Eq. ^ with 
boundary conditions tp{±L,t) = po for large values of L. 
Typical outcomes are summarized in Fig. [TJ Comparing 
the evolution of the density at a; = for different widths 
of the defect [Fig.[lja)] we find that, for sufficiently small 
^, n(0, t) decays monotonically, while non monotonic be- 
havior appears for C. > Ccr ~ 2.26. The dynamics allways 
ends up in some asymptotic regime, and the differences 
in the ways of how this regime is achieved are due to the 
excitation of sound waves, which are clearly seen in the 
case of the narrow defect [Fig. [He)], while strongly atten- 
uated for the case of the wide defect [Fig. [T{d)] . We have 
also simulated the evolution of an initially uniform BEC 
when the defect is conservative [Fig. [Tfe)] While 
the creation of a minimum of the atom density and the 
excitation of sound waves propagating outwards the de- 
fect look similar to the case of the dissipative defect, the 
current density reveals significant differences [Fig. [TI[f)]. 
In the former case, the currents are transient and due to 
the kicking of the atoms from the high-potential regions 
around the defect. However, the dissipative defect is a 
sink of energy thus, the currents point towards the defect 
and are permanent in this scenario. The presence of the 
dissipative defect implies the imprinting of a nontrivial 
phase on the condensate wavefunction [Fig. [iff)]. 

Fig. Wijo) shows another interesting phenomenon: a 
more intensive defect results in less intensive loss of con- 
densed atoms at large times. More efficient atomic evap- 
oration by a weaker defect can be explained from the 



fact that such defect created a hole with a finite density 
leading to the long-term action of the dissipation, while 
a strong defect results in fast local density decay to zero 
and thus to suppressed long-term effect of the dissipation. 




FIG. 1: (a) Evolution of the density for different widths of 
the defect ((3]) with 70 = 0.5. The dashed red line shows the 
exponential decay (b) Dynamics of N{t) normalized to its 
initial value A^o = 2poL for different strengths of the defect. 
In (c) and (d) surface plots of n for the defect with = 1 < ("cr 
and = 5 > ^cr , respectively. In (e) the same as in (d) but for 
70 — ~0.5i (conservative defect), (f) Snapshots of the current 
density j — — IV'j^ • [arg(V')]i at t — 100 for the defect from 
(d) (black solid line) and for the conservative defect from (e) 
(red dashed line). In all panels a — 1 and po — 0.4. 

Generation of dark solitons.- Next, we discuss the en- 
suing dynamics after removing the dissipative defect at 
time tf (Fig. [5]). It can be seen how the "initial" dark 
hole splits into a packet of sound waves and dark soli- 
tons propagating outwards the impurity. This part of 
the dynamics is described by the integrable NLS equa- 
tion allowing for exact computation of the number of 
solitons and in this sense is not related to the type of 
the defect creating the initial intensity hole. The main 
difference between conservative and dissipative defects, 
however, steam form the initial phase shown in panels (c) 
and (e) of Fig. [TJ in the last case the excitation of linear 
waves (and eventually small amplitude dark solitons) is 
strongly suppressed, what makes the process more con- 
trollable. A qualitative feature of the phenomenon dis- 
played in Fig. [1] (f) is that while the dissipative defect 
is acting, its main effect is imprinting a phase variation 
propagating outwards the defect. In Fig. [2] one observes 
the small amplitude waves generated at i < as dark 
solitons created by the phase imprinting, while at when 
the defect is switched off, the shape of the dark hole 
changes generating two well-separated deep (slow) dark 
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solitons (at this stage, the generation of sohtons is dom- 
inated by the density engineering). In other words the 
dissipative defect provides a controllable way of generat- 
ing dark solitons based on phase and density engineering. 



(a) n{x,t) |k (b) ^ yc) 



FIG. 2: Generation of dark solitons for 70 = 0.5, po ~ 0.4 
and a — 1. Shown are density plots with (a) = 0.1, = 15; 
(b) C = 0.5, = 15 and (c) ( = 0.5,*/ = 40. The time range 
shown is t £ [0, 100] and the spatial range x G [—50, 50]. 



The situation in a trap might be very different be- 
cause of a finite number of condensed atoms, resulting 
in degradation of the whole picture after significant loss 
of particles. To study this scenario, we have considered 
an additional term in Eq. ([T]) of the form v'^x'^tp where ly 
determines the linear oscillator frequency of the longitu- 
dinal harmonic trap. In Fig. EI a) we present the evolu- 
tion of the density in the center of the trap for different 
widths of the dissipative defect. The initial dynamics 
(not visible on the scale of the figure) is very similar to 
that of the case with constant background [Fig. [BJa)]. 
With time the loss of particles results in a significant de- 
crease of the background what takes the system beyond 
the mean-field regime. This happens however after a suf- 
ficiently long time: « 800 in our dimensionless units. 
When the dissipative effect is switched off soon enough. 
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FIG. 3: Effect of the dissipative defect on the condensate in a 
trap ior a = 1, u ^ 0.002, po = 0.4, and 70 = 0.5. (a) n(0, t) 
for different defect widths: = 0.5; 2; 5. (b) Pseudocolor plot 
of the density evolution for — 2. (c) Same as in (b) but 
removing dissipation ai tf — 100. 

we still observe [Fig.[3](c)] the generation of dark solitons 
which are later reflected by the parabolic potential and 
start their characteristic oscillatory motion 0, [l^ • 

Bright solitons.- When the interatomic interactions are 
attractive, cr < 0, a uniform background is modulation- 
ally unstable. In that case, the instability is triggered 
by the dissipative defect leading to the emergence of the 



spatial structures as it is shown in Fig. [4] In compari- 
son with the repulsive case we observe a faster decay of 
the amplitude in the defect domain due to the stronger 
attraction of atoms by the populated regions than by 
the defect region (in the repulsive case atoms are re- 
pelled from the densely populated domains). The devel- 
opment of the instability leads to the creation of solitonic 
pulses [Fig.[4](c)]. This is another counter-intuitive phe- 
nomenon since solitons are characteristic excitations of 
Hamiltonian conservative systems. However, it is easy to 
see that after the density hole is created the loss of the 
atoms is significantly reduced and almost everywhere the 
system remains to be "almost-conservative" . 




FIG. 4: Generation of bright solitons for 70 = 0.5, po = 0.4. 
(a) n(0, i) for different defect widths — 1,3,5). (b) density 
profile at i = 50 (b) and (c) density evolution for C = 3- 
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FIG. 5: Generation of dark (a) and bright (b) gap solitons 
after switching off the defect, with 70 — 0.5 and = 2, at 
tf = 40 and t/ = 100, respectively (solid lines). The initial 
profiles (dashed lines) are Bloch states for A = —1 and po = 
0.5. The dashed line indicates the defect profile. 

Gap solitons.- Optical lattices (OLs) have proven to 
be one of the most successful tools for matter-wave man- 
agement [l^]- Now we consider again Eq. dJ) but with 
the addition of the periodic potential V{x) = Acos{2x), 
whose period is chosen to be tt provided the amplitude A 
is measured in the units of the recoil energy Er. The 
background solution must be taken as a Bloch state, 
ipa (x) bordering an edge of a gap of the energy spectrum, 
'ip{x,0) = po^Paix) {po measures the density amplitude 
and we consider the first lowest gap). 

We consider a — —1 and moderate OL depths {A ~ 1). 
Then a Bloch state background is stable (unstable) when 
it borders the bottom (top) of a band corresponding to 
positive (negative) effective mass. The behavior of Bloch 
states affected by the localized dissipative defect closely 
resemble the dynamics in a homogeneous condensate, 
with the major qualitative difference that both dark and 
bright solitons can be generated in an attractive BEG, 
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FIG. 6: Initial (thin) and final (thick) at f = 1000 profiles 
of the density with, a = —1, (a) and without, a = 0, (b) 
the nonlinearity. The parameters are po = 0.4, ^ — 0.5, and 
A — —10. The dashed line illustrates the defect. 

starting with states with positive and negative effective 
masses, respectively [Fig. [5]. 

In deep optical lattices |^| > 10, however, we have ar- 
rays of well separated condensates, among which atomic 
exchange can occur due to the tunneling and the tight 
binding approximation becomes more adequate [l8| . The 
most exciting effect which can be observed in this situ- 
ation is the creation of persisting density holes due to 
the suppression of tunneling by the attractive nonlinear- 
ity. Fig. [6l where we compare the dynamics with and 
without the nonlinearity, illustrates this phenomenon. 




FIG. 7: Density plots on the region {x,y) £ [-30,30] x 
[—30,30] for different times. The parameters are tf — 10, 
70 = 0.5, C = 2, po = 0.66, and u = 0.01. 

The stationary density hole created in the case cr = — 1 
(panel a) is not a dark gap soliton: in our case the initial 
alternating phase of the Bloch state is preserved and the 
concept of effective mass loses its significance because 
the bands are very narrow. Thus the created hole is 
supported only by the interplay between the tunneling 
and the attractive nonlinearity and does not have any 
analogy in the case of a homogeneous medium. The phe- 
nomenon can be understood in terms of the effective po- 
tential. Indeed, after long enough time the distribution 
n{x) becomes (quasi-)stationary, even in the presence of 
the dissipative defect, and Vcg{x) = V{x) — n{x) can 
be interpreted as the effective potential affecting a sin- 
gle particle. Since for V ^ 10 the band width is of 
order of O.QlEr and n ^ 0.15 Er, in the absolute mini- 
mum of Vcs the lowest energy levels are shifted down, the 
shift being an order of magnitude larger than the energy 
splitting due to hopping, resulting in the suppression of 
tunneling to the central cell. 

Multidimensional examples.- A localized dissipative 
defect can be used for generating higher dimensional non- 



linear structures. As an example we have simulated the 
effect of a 2D Gaussian defect ([3]) in a repulsive pancake 
shape BEG in a harmonic trap i^^{x^ +y'^). After switch- 
ing off the defect at time tf we observe the generation of 
ring dark solitons [l^ , as shown in Fig. [T] 

Conclusion.- We have shown how a dissipation source 
can be used to generate nonlinear coherent excitations 
such as dark, bright, gap or ring dark solitons by en- 
gineering both the phase and amplitude in systems de- 
scribed by NLS equations. Our results can be apphed to 
matter wave management in BEGs. 

V.A.B. acknowledges support from the FGT grant, 
PTDG/FIS/64647/2006. V.M.P.-G. is supported by 
grants FIS2006-04190 (Ministerio de Giencia e Inno- 
vacion, Spain) and PGI-08-093 (Junta de Comunidades 
de Gastilla-La Mancha). 



[1] S. Burger et ai, Phys. Rev. Lett. 83, 5198 (1999). 

[2] S. Stellmer et ai, Phys. Rev. Lett. 101, 120406 (2008); 

A. Weller, et al, Phys. Rev. Lett. 101, 130401 (2008). 
[3] K. E. Strecker, G. B. Partridge, A .G. Truscott, and R. 

G. Hulet, Nature 417, 150 (2002); L. Khaykovich et al., 

Science 296, 1290 (2002). 
[4] S. L. Cornish, S. T. Thompson, and C. E. Wieman, Phys. 

Rev. Lett. 96, 170401 (2006). 
[5] M.R. Matthews, et ai, Phys. Rev. Lett. 83, 2498 (1999); 

J.R. Abo-Shaeer, et al, Science 292, 476 (2001); V. 

Schweikhard, et al, Phys. Rev. Lett. 93, 210403 (2004). 
[6] B. P. Anderson et al, Phys. Rev. Lett. 86 2926 (2001). 
[7] B. Eiermann et ai, Phys. Rev. Lett. 92, 230401 (2004). 
[8] Z. Dutton et al, Science, 293, 663 (2001); J.J. Chang, P. 

Engels, and M. A. Hoefer, Phys. Rev. Lett. 101, 170404 

(2008). 

[9] N.S. Ginsberg, J. Brand and L.V. Han, Phys. Rev. Lett. 
94, 040403 (2005); C. Becker et ai, Nature Phys. 4, 496 
(2008). 

[10] P. Engels, C. Atherton, and M. A. Hoefer, Phys. Rev. 

Lett. 98, 095301 (2007). 
[11] E. A. Burt, et al, Phys. Rev. Lett. 79, 337 (1997). 
[12] M.-O. Mewes, et ai, Phys. Rev. Lett. 77, 988 (1996); 

J.R. Abo-Shaeer, C. Raman, and W. Ketterle, Phys. Rev. 

Lett. 88, 070409 (2002). 
[13] T. Lahaye et al, Phys. Rev. Lett. 101, 080401 (2008). 
[14] A. Alexandrescu, V. M. Perez-Garci'a, Phys. Rev. A 73, 

053610 (2006); R. Pardo, V. M. Perez-Garcfa, Phys. Rev. 

Lett. 97, 254101 (2006). 
[15] T. Gericke, et ai, Nature Phys. 4, 949 (2008). 
[16] V. A. Brazhnyi and V. V. Konotop, Phys. Rev. A 68, 

043613 (2003). 

[17] The data in panels (e), (f) correspond to the experiment 

by M. R. Andrews, et al. Phys. Rev. Lett. 79, 553 (1997). 
[18] see e.g. V.A. Brazhnyi and V.V. Konotop, Mod. Phys. 

Lett. B 18, 627 (2004); O. Morsch and M. Oberthaler, 

Rev. Mod. Phys. 78, 179 (2006). 
[19] Y. S. Kivshar, X. P. Yang, Phys. Rev. E 50, R40 (1994); 

P. G. Kevrekidis, et al, Phys. Rev. Lett. 90, 230401 

(2003). 



